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ABSTRACT 

The last few years have seen a variety of new image mask designs for 
diffraction-hmited coronagraphy. Could there still be useful designs as yet undis- 
covered? To begin to answer this question, I survey and unify the Fraunhofer 
theory of coronagraph image masks in the context of a one- dimensional classical 
coronagraph. I display a complete solution to the problem of removing on-axis 
light assuming an unapodized entrance aperture and I introduce the attenua- 
tion function, a measure of a generic coronagraph's off-axis performance. With 
these tools, I demonstrate that the masks proposed so far form a nearly complete 
library of image masks that are useful for detecting faint extrasolar planets. 

Subject headings: astrobiology — circumstellar matter — instrumentation: adap- 
tive optics — planetary systems 

1. INTRODUCTION 

In a classical coronagraph (Lyot 1939), an image mask placed at a telescope's focus and 
a Lyot stop in a succeeding pupil plane block most of the light from a bright on-axis source 
allowing the telescope to better image faint off-axis sources. Classical coronagraphs have 
provided the first images of a brown dwarf orbiting a nearby star (Nakajima et al. 1995) 
and striking images of circumstellar disks (e.g. Smith & Terrile 1984; Clampin et al. 2003) 
and promise to enable imaging of bright extrasolar planets within the next decade (Trauger 
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et al. 2003). Initial studies suggest that classical coronagraphs can potentially image even 
extrasolar terrestrial planets (see the review by Kuchner & Spergel 2003). 

The last few years have seen a cornucopia of new coronagraph designs, including a 
variety of new image mask designs, typically invented using Praunhofer diffraction theory. 
Perhaps construction tolerances and vector electromagnetic effects will limit the practicality 
of all these designs. However, the more choices of masks we have explored using Fraun- 
hofer theory, the more likely it seems we may find a useful one for terrestrial planet finding. 
With extrasolar terrestrial planets in mind, I attempt to develop a unified picture of corona- 
graph image masks to determine the range of potentially useful solutions to the Fraunhofer 
diffracted-light problem. I concentrate first on a one-dimensional coronagraph with a tophat 
entrance aperture and then consider more general entrance apertures. 



We will examine a one-dimensional coronagraph (Sivaramakrishnan et al. 2001) com- 
prising an entrance aperture. A, an image mask, M, and a Lyot stop, L, each of which is 
represented by a complex- valued function constrained to have absolute value < 1. We use 
the notational conventions of Kuchner & Traub (2002) and Kuchner & Spergel (2003); letters 
with hats represent pupil plane quantities, and quantities transform as follows: 



We will assume the optics have removed any quadratic phase terms associated with small 
focal lengths. We will set all time- varying factors in the electric fields equal to 1; whenever 
we say "field", we mean the amplitude of the time varying field. Mostly we will discuss 
monochromatic effects, though sometimes we will consider the effects of finite bandwidth. 

We will use dimensionless units, but our units can be translated into physical distances 
as follows. In the pupil planes, our coordinates {u, ui) measure distance from the optical axis 
in units of the local pupil diameter, which for example, at the primary mirror would simply 
be D, the diameter of the primary mirror. However, since the model is one dimensional, 
the primary mirror doesn't have a diameter so much as a width; it extends forever in one 
direction. In the image planes, our coordinates {x, 9, etc.) measure distance from the optical 
axis in units of the diffraction scale, which is ordinarily A/, where A is the wavelength of 
light and / is the focal ratio. In the plane of the sky, they measure angle from the optical 
axis in units oi X/ D radians. 

An incoming wave incident on the entrance aperture creates a field with amplitude 
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E{u). In our model, when an incoming wave interacts with a stop or mask, the function 
representing the mask multiphes the wave's complex amplitude. So after the wave interacts 
with the entrance aperture, the amplitude becomes A{u) ■ E{u). 

After interacting with the entrance aperture, the beam propagates to an image plane, 
where the new field amplitude is the Fourier transform of the pupil plane field amplitude, 

A{x) * E{x), where * denotes convolution. There the beam interacts with the image mask, 
and the field amplitude becomes M{x) ■ {A{x) * E{x)). Then the beam propagates back to 
a second pupil plane, where the field amplitude is M{u) * {A{u) ■ E{u)). 

In the second pupil plane, the wave interacts with a Lyot stop, changing the field 
amplitude to L{u) ■ [M{u) * {A{u) ■ E{u))]. Then the beam propagates to a final image plane, 
where the final image is L{x) * [M{x) ■ {A{x) * E{x))]. For a point source, the intensity of 
the final image is proportional to the absolute value of this quantity squared. 

A coronagraph aims to reduce the final image of an on-axis point source, for which E(u) 
is a constant, which we can set equal to 1. The field created by such a source in the second 
pupil plane, after the Lyot stop, is L{u) ■ (M * A) (Sivaramakrishnan et al. 2001). The 
corresponding final image field is the Fourier transform of this quantity, FT(L • (M * A)). 
Masks for which 

[L ■ {M * A)]{u) = (2) 

will block identically all of the light from an on-axis point source; we aim to find these masks. 
This paper concerns solutions and approximate solutions to this linear problem. 



2.1. Tophat Entrance Aperture, Arbitrciry Lyot Stop 

We will examine first a system where the entrance aperture, A{u), is opaque {A{u) = 0) 
for \u\ > 1 and transparent {A{u) = 1) for < 1/2, and the Lyot stop, L{u), is opaque 
for \u\ > (1 — e)/2, where e < 1. To understand Equation 2 in the context of our simple 
coronagraph, we can write A{u) as a difference of two Heaviside step functions, H{u): 

A{u)^ H{u + 1/2)- H{u- 1/2) (3) 

Then, since convolution with a Heaviside step function represents indefinite integration, we 
can write 

M{u) * A{u) = M{u+ 1/2) -M{u- 1/2) (4) 
where {d/ du)M.{u) = M{u). Now we can see that Equation 2 demands only that 

M{u + 1 /2) ^ M{u- 1/2) for -(1 - e)/2 < < (1 - e)/2 (5) 



-4- 



or equivalently. 

M{u)^M{u-l) for e/2 < < 1 - e/2. (6) 

An uncountable infinity of non-trivial masks meet this symmetry criterion. Figure la 
shows a generic function, A4{u) that satisfies Equation 6. In the region e/2 < u < 1 — e/2, 
M{u) looks the same as it does in the region —e/2 < u < —(1 — e/2); outside those regions, 
the function shown in Figure la has some random squiggles, to remind you that Equation 6 
doesn't specify anything about the function there. 



2.2. Notch Filter Masks and Band-Limited Masks 

Figure lb shows the simplest interesting solution to Equation 6: 

At (li) = constant for e/2 < < 1 - e/2. (7) 
This statement about M.{u) translates into two requirements on M{u): 

M{u) = for e/2 < \u\ < 1 - e/2 (8) 

and 

l-e/2 

/ M{u) du = 0. (9) 

J-e/2 

We can call functions N{x) whose Fourier transforms, N{u), satisfy Equation 8 "notch filter" 
functions. Masks satisfying Equation 7 or equivalently, both Equations 8 and 9, are called 
notch filter masks (Kuchner & Spergel 2003; Debes et al. 2004). 

Notch filter masks for which M{u) = for \u\ > e/2 are called band-limited masks 
(Kuchner & Traub 2002). Band-limited masks are an important subset of notch- filter masks 
because as we will see below, only the low-wavenumber components of any mask affect the 
off-axis light. Naturally, band-limited functions display wide variation, from the M{x) = 
sm^{nex/2) mask that optimizes inner working angle, to masks where M{x) is a prolate 
spheroidal wavefunction (Kasdin et al. 2003). These latter masks optimize search area. 
The Gaussian mask (Sivaramakrishnan et al. 2001) is a good approximation to the prolate- 
spheroidal wavefunction band-limited mask. 




Fig. 1. — Ai{u) for a) a generic solution of Equation 6, b) a zeroth order, or notch filter 
mask c) a first order mask d) a 1-D disk phase mask, a good approximation to a first order 
mask. Eliminating on-axis light in a 1-D coronagraph with tophat entrance aperture requires 
finding a function with this translational symmetry. 
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2.3. First Order Masks 

The next simplest solution to Equation 6 is one that adds a linear slope to A4{u), while 
maintaining the needed translational symmetry. I.e., 

.., . ( kiu for e/2 < M < 1 - e/2 

■^("'"U.l^+D £or-(l-e/2)<«<-e/2. ('"^ 
Figure Ic shows an example of such a sawtooth-hke function. When M{u) takes this form, 

M{u) = ki for e/2 <\u\<l- e/2 (11) 

and 

/ M{u)du= -ki{l-e). (12) 

J-e/2 

Equations 11 and 12 contain Equations 8 and 9 as a special case (fci = 0). 

Mask functions M{x) that solve Equations 11 and 12 can always be decomposed into a 
sum of a notch-filter function and a 5-function. For example, a 1 — 5-function mask 

M{x) = l-5{x) (13) 
M{u) = 5{u)-l (14) 
M{u) = H{u)-u (15) 

fills the bill. We can call notch filter masks "zeroth order" masks, and this 1 — 5-function 
mask a "first order" mask. 

Of course, we can not biiild this 1 — 5— function mask. The closest we can come is to 
build a mask that is everywhere equally transparent, but that generates a phase shift of tt 
in the center: 

N f 1 for \x\ > 1/4 

^ ' [ -1 for < 1/4. ^ ' 

M{u) = 6{u) -sm{7ru / 2)/ {nu/ 2) (17) 

M{u) ^ H{u)-u + —u^ - ... (18) 

Figure Id shows Ai{u) for this mask — it approximates the sawtooth-like exact solution of 
Figure Ic. The disk phase mask of Roddier & Roddier (1997) is a two-dimensional version 
of this approximate solution. 
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2.4. The Rest of the Series 

Equation 6 shows that over some interval, Ai{u), and therefore M[u) are periodic 
functions, with period 1. We can expand such functions as a Fourier series: 

oo 

P[u) = ^ Aj sin 271 jx + Bj cos 27rjx (19) 

j=0 

where Aj and Bj are complex. This series satisfies Equation 6. However, this solution would 
completely specify the function over all x; Equation 6 offers more freedom than that. We 
can express a complete solution as the sum of the above Fourier series and any notch filter 
function: 

oo 

M{u) = N{u) + Aj cos 27r(j - l)x + Bj sin 27r(j - l)x (20) 
i=i 

oo 

M{x) = N{x) + J2CjS{x-{j-l)) + DjS{x + {j-l)), (21) 

i=i 

where and Cj — {Aj — iBj)/2 and Dj — {Aj + iBj)/2, as long as 

I M{u)du = / M{u)du (22) 

J-e/2 Je/2 

— £^4^ sinc7r(j — l)e. (23) 

In other words, all masks that completely remove on-axis light in a coronagraph with a 
tophat entrance aperture can be represented as the sum of a notch filter function and a 

series of 5-functions located x — . . . , —2, —1, 0, 1,2, So far, we have only considered 

masks with all Aj and Bj equal to 0, which we called zeroth order masks, or masks with 
only Ai 7^ 0, which we called first order masks (there is no B-i). We can refer to masks with 
Aj or Bj ^ 0, for j > 1 as "higher order" masks. 

We can also consider the situation where the Lyot stop is bigger than the entrance 
aperture; e < 0. In this case, the solution to Equation 2 takes the form 

oo 

M{u) = H{u) + Aj cos 2-Kjx + Bj sin 27: jx (24) 

oo 

M{x) = H{:,) + Y,Cj5{x-{3-l)) + Dj5{x+{3-l)), (25) 

where H{x) is a high-pass filter function, i.e. one where H{u) — for |ti| < 1 — e. High-pass 
filter functions do not have any effect on the image, however; the functional part of the mask 



-8- 



is the series of 5-functions. One example of a mask that works (approximately) for e < is 
the 1-D disk phase mask illustrated above. It does not suddenly fail if e becomes small or 
negative — though the approximation to a (5-f unction degrades towards the edges of the pupil 
plane. 

Masks containing a row of (5-functions may not be buildable, but we could build approxi- 
mate versions, like the disk phase mask, and we could build low-pass-filtered versions of some 
masks by convolving the mask function with sinc(27ra;) for example, (though constructing 
such masks could be tricky!). 

However, all masks in this series except for notch filter masks and first order masks face 
the following complication in a broad-band system. In a real optical system, M{x) necessarily 
changes with the wavelength of the light. A; for an intensity mask, M{x) ^ M{xX/\o), where 
Aq is the reference wavelength used for the purpose of defining the diffraction scale. Other, 
wavelength-dependent effects may appear, or be built into the mask, but this dilation is 
intrinsic to a classical coronagraph. 

Consequently, Ai{u) ~ Ai{uXo/ X). The notch-filter solution, for which Ai{u) is con- 
stant over a large range in u, can trivially be made to work over a large range in A because of 
its dilation symmetry. However, all other solutions will only work at one wavelength, unless 
M{x) explicitly changes as a function of wavelength in some special way other than simple 
dilation. 



3. OFF-AXIS LIGHT AND A MORE GENERAL CORONAGRAPH 

Now we will assess the impact of a coronagraph on an off-axis source, like an cxtrasolar 
planets located at angle 6 from the optical axis in the plane of the sky. We can proceed in 
the context of a more general problem, where A{u) and L{u) can be any functions at all that 
we can Fourier transform. 



3.1. The Attenuation Function: Inner Working Angle 

A band-limited mask is what-you-scc-is-what-you-gct; a coronagraph with a band- 
limited mask attenuates the field from a distant source by a factor of M(9) and the Point 
Spread Function (PSF) in such a coronagraph does not depend on 9. (Kuchner & Traub 
2002). But in general, the PSF in a coronagraph depends on 9, and clearly, the shapes of 
many image masks do not directly show how the mask attenuates off-axis sources. How can 
we understand in general what a coronagraph does to image of an astronomical source? 
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Kuchner & Traub (2002) showed that the high wavenumber components of the mask 
function, those with > 1 — e/2, do not affect how the mask interacts with off-axis hght. 
Consequently, for any mask, we can convolve M{x) with sinc(27r(l — e/2)x) and get the same 
effective mask. For a notch filter mask, this operation just means looking the band-limited 
part of the mask, which as we mentioned, is trivial to interpret. We will apply a similar idea 
here to reveal the workings of any image mask. 

Consider a point source, providing a field 6{x—Xi) in the plane of the sky, and a harmonic 
mask function M{u) = m5{u — Ui). The field after the entrance pupil is A{u) exp(— 27rmxi), 
and the field in the first image plane is A{x — Xi). The field after the image mask is 
mey.]i{2'Kiuix)A{x — xi). The field in the second pupil plane is mA{u — Ui) exp[—2ni{u — 
Ui)xi) The field after the Lyot stop is mL{u)A{u — Ui) exp(— 27ri(M — Ui)xi). The final field 
is mex^{2muiX-i)¥T[L{u)A{u — u\)\ *5{x — xi). In other words, the point source response of 
this particular special coronagraph is an intensity pattern \m¥T[L{u)A{u — ui)\^ centered 
&X, x — xi. 

We can evaluate the field amplitude at the center of the pattern simply by replacing 
the Fourier transform with Jf^ du. Then we see that the amplitude attenuation provided 
by this coronagraph at the center of the image of a point source is 

F{xi) = me^""'^i [L{ui) * A{-ui)] for a harmonic mask. (26) 

The intensity at the center of the Point Spread Function (PSF) is \F{xi)\^. We can call F{xi) 
the coronagraph's amplitude attenuation function and |F(xi)|^ the coronagraph's intensity 
attenuation function. 

We can find the attenuation function for any arbitrary mask by expanding the mask 
function as M{u) — Jf^ M{ui)5{u — ui)du. Then the attenuation function becomes a 
Fourier integral: 

/oo 
M{u) e^™^^ [L{u) * A{-u)] du (27) 
-oo 

= ¥T{M ■ {L{u) * A{-u))){xi) (28) 
= M*[l-A'^^ (xi), (29) 

where f denotes complex conjugation. This definition applies for any arbitrary entrance 
aperture, image mask, or Lyot stop. The PSF shape may change as a function of Xi, the 
location of the off-axis source. However, the center of the PSF is usually the maximum of 
the PSF, so the attenuation function measures a critical feature of the mask for the purpose 
of finding off-axis point sources, like extrasolar planets. 
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Equation 29 shows that while a mask may look to us like it has sharp edges, in some 
sense it looks to the sky like it is smooth — convolved with a smoothing kernal of the form 
L ■ . For example, the 5— function mask described above has an amplitude attenuation 
function F[x) = 1 — [L ■ A^lx). If the Lyot stop is a tophat function and e is small, then 
F{x) 1 — sinc^(7ra;/2). If the image mask is band-limited, then F{xi) = M{xi). 

Figure 2 shows how F{x) is constructed for a coronagraph with a tophat entrance 
aperture, tophat image mask, and tophat Lyot stop. The left column shows image plane 
quantities, and the right column shows their pupil plane conjugates. Figure 2a, b, and c 
show the entrance aperture, A. the image mask, M, and the Lyot stop, L. Figure 2d shows 
the convolution of L{u) and A{u) and the image plane conjugate of this function, L{x)A{x). 
These quantities act respectively as filter function and smoothing kernal. The amplitude 
attenuation function, F{x), shown in Figure 2e, is a smoothed version of the mask. 

In this example, L{u) * A{u) is a constant for \u\ < e/2, and zero for |m| > 1 — e/2, as 
shown in Figure 2d. Consequently, if M{u) were a band-limited mask function with power 
only where \u\ < e/2, then multiplying M{u) by L{u) * A{u) would only serve to multiply 
M{u) by a constant. In other words, M{x) would be an eigenfunction of convolution with 
L(x)-A(u). 

Figure 3 shows M{x) and \og\F{x)\'^ for several different image masks: a phase knife 
(Equation 31), a tophat mask, a one-dimensional disk phase mask (Equation 17), a Gaussian 
mask, and a band-hmited mask (M{x) = 1 — sinc^(7rex/2)). Of the masks shown, only the 
phase knife and the band-hmited mask provide F(0) = 0. Of those, only the band-hmited 
mask provides perfect cancellation of on-axis light and the deep null needed for terrestrial 
planet finding. 

The attenuation function provides a good metric for the inner working angle, djw of 
a coronagraph. For example, we could define the inner working angle by \F{9iw)\'^ = 1/2. 
With this definition, and assuming a tophat Lyot stop with e = 0, 9iw = 0.58\/D for 
the 5-function image mask, and 0A3X/D for a phase knife. For comparison, the smallest 
possible inner working angle for a band-limited mask is 6jw = 0.6AX/D — except that for 
this band-limited mask, e = 1, so the coronagraph would have no throughput. 

Notice that with a mask that is not a notch filter mask, apodizing the Lyot stop makes 
L into a broader function, which increases the coronagraph's inner working angle. Likewise, 
we could push the inner working angle inwards if we replaced the Lyot stop with a pair of 
pinholes ai u = ±1/2, for example. However, in reality, these distinctions probably don't 
matter, because low order aberrations, like pointing error, will probably set the inner working 
angle of a real coronagraph 
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Fig. 2. — Constructing the attenuation function, a smoothed version of the mask that the 
planets see. a) entrance pupil and its Fourier conjugate b) image mask and its conjugate c) 
Lyot stop and its conjugate d) the smoothing kernal and the filter function e) the amplitude 
attenuation function, F{x), and its conjugate. 




Fig. 3. — Mask functions, M{x), and intensity attenuation functions, |F(a;)p, for several 
image masks, calculated assuming the entrance aperture and Lyot stop are identical tophat 
functions, a) Phase Knife b) Tophat c) 1-D Phase Disk d) Gaussian e) 1 — sinc^ (Band- 
Limited) 
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3.2. Mask Symmetry and Stellar Leak 

By taking derivatives of Equation 29, we can construct a Taylor expansion for F{x) 
about X = 0: 

— / Miu) i2muxiY [Liu) * A(-u)] du. (30) 

Consider a system where L{u) and A{u) have even symmetry. If M{u) has even symmetry 
then (^) " F{x) will vanish for all odd n. If M{u) has odd symmetry then (^) " F{x) 

will vanish for all even n. For example, in a coronagraph where A{x) and L{x) have even 
symmetry, if we use any image mask with odd symmetry, the on-axis final image field from 
an on-axis point source will be zero. 

Mask functions, M(a;), with odd symmetry necessarily become negative, so they require 
manipulating the phase of the beam. The simplest such mask is the phase-knife coronagraph 
(Abe et al. 2001), for which 

1 for a; > 

M{x) = { Q for X = (31) 
—1 for a; < 0. 

A phase knife consists of a half-plane of glass joined to a half-plane of phase-retarding 
material so that the seam falls on the optical axis. We may be able to make a more practical 
version of this mask by laying an opaque strip over the seam so that the mask function 
becomes 

{1 for x > s 
for \x\ < s (32) 
—1 for X < —s. 

where s is one or two diffraction widths. This variation retains the mask's odd symmetry, 
but removes the seams from direct illumination. 

Though F{x = 0) = for the phase knife, the mask does not ehminate all on-axis light. 
The on-axis source may produce an image with zero central intensity, but with bright wings 
that can overwhelm the image of an off- axis planet. A perpendicular pair of phase knives, 
hke the four-quadrant phase mask (Rouan et al. 2000; Riaud et al. 2001, 2003; Lloyd et al. 
2003) does a better job than a single phase knife at reducing these wings at the cost of some 
search area and also exactly satisfies Equation 2 for a circular entrance aperture and circular 
Lyot stop. 

However, the wings of the image of an on-axis point source are not the only drawback 
to the phase-knife mask. Real stars do not provide perfect on-axis point sources; stars have 
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Tablc 1: Null Depth 
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"Assuming appropriate pupil stops that the zeroth order leak is negligible. 

finite angular size and telescopes can not point at them perfectly accurately (e.g. Riaud et al. 
2001; Kuchner & Traub 2002; Kuchner & Spergcl 2003; Lloyd et al. 2003). We must consider 
the intensity leak from slightly off-axis starlight, which we can approximate by \F{9)\'^, where 
9 is the angle off axis. 

A common misconception is that disk phase masks are especially sensitive to pointing 
error. However, any coronagraph has a pointing error tolerance related to its inner working 
angle; the smaller the inner working angle, the tighter the tolerance. Disk plase masks have 
small inner working angles, so they have tight pointing tolerances. But at a given inner 
working angle, disk phase masks perform well. 

In general, the leak increases as some power of 9. Equation 30 shows that the intensity 
leak from a mask with odd symmetry, like a phase-knife mask, is 0{9'^). A notch filter mask 
has even symmetry, so the intensity leak from a notch-filter mask is 0{9'^). The circular 
phase mask (Roddier & Roddier 1997) and the Gaussian mask have F{x = 0) 7^ 0, so they 
leak at 0{9^), though the zeroth order leak can be made neghgible. All masks of containing 
sin terms in Equation 21 produce an 0(9"^) intensity leak. 

Viewed at quadrature, reflected visible light from the Earth is 2 x 10~^° times as bright 
as direct hght the Sun. We do not need to suppress the starhght to quite this level in the 
center of the image plane since the planet appears typically a few diffraction widths from 
the star in the wings of the stellar leakage, which we can control with an apodized Lyot stop 
if necessary. A reasonable assumption might be that for terrestrial planet finding, we can 
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tolerate up to F{A9) — 10 ^, where is the effective pointing error (which may contain 
some contribution from the star's finite angular size). 

Given this assumption, we must center the star on the mask to an accuracy of roughly 
= 10-^/ i^Oiw, where the mask's leak is and Oiw is the coronagraph's inner work- 

ing angle. For example, a fourth-order intensity leak, or a fourth-order "null", to use the 
language of interferomctry, translates into a pointing requirement of A^ < ^/i^/lOO. Table 1 
summarizes the leading order intensity leaks from the masks discussed in this paper. 

3.3. Other Directions: Apodized Entrance Apertures 

We found a complete solution to the problem of removing on-axis light in a one- 
dimensional coronagraph with an un-apodized entrance aperture. However, several coro- 
nagraph designs, classical and otherwise, use apodized entrance pupils (Kasdin et al. 2003) 
or specially shaped pupils (Spergel 2001; Kasdin et al. 2003; Vanderbei et al. 2003a,b). 
Other designs use a pair of shaped mirrors to generate an apodized beam (Guyon 2003; 
Traub & Vanderbei 2003). These designs aim to suppress the wings of the PSF so far that 
the edges of the image mask are not illuminated by the image of an on-axis source, so a 
simple tophat mask suffices to remove the on-axis fight to the necessary level. Achieving 
this suppression generally requires numerical optimization to find an approximate solution 
to Equation 2. Goronagraphs with apodized entrance apertures, like coronagraphs with 
apodized Lyot stops, are especially robust to low-order aberrations. 

We can also find new exact solutions to Equation 2 using an apodized entrance aperture 
and a phase mask (Aime et al. 2002; Soummer et al. 2003). However, these exact solutions 
require M{x) to vary in a complicated way with wavelength. One design that provides a 
self-adjusting M(x, A) to be combined with an apodized entrance aperture is the achromatic 
dual-zone phase coronagraph (Soummer et al. 2003). Achromatic dual-zone phase masks can 
be designed with (9(^^^) leak. Combining image masks and apodized entrance pupils remains 
relatively unexplored territory — though perhaps such designs shold not be called classical 
coronagraphs. 

4. CONCLUSION 

According to Fraunhofer diffraction theory, uncountably many non-trivial image mask 
designs can completely remove the light from an on-axis source in a classical coronagraph 
with an un-apodized entrance aperture. We organize these designs using a trigonometric 
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series solution to the equation L ■ [M * A) — 0. Notch filter masks are the zeroth order 
solutions, and the only solutions that are trivially achromatic in our approximation over a 
large bandwidth. Higher order solutions can be decomposed into the sum of a notch filter 
function and a series of 5{x) functions; the the disk phase mask (Roddier & Roddier 1997) 
approximates one term in the sum. 

We defined the mask amplitude attenuation function, F{x), as the amplitude at center 
of the final image of a point source located at angle x off-axis in the plane of the sky. This 
function, which we showed is a smoothed version of the mask function, provides an easy 
way to judge the null depth and inner working angle for non-band-limited masks. Through 
this analysis, we showed that only masks with even symmetry can provide the null depth 
necessary for terrestrial planet searches; for example, masks based on a phase knife do not 
provide the necessary null depth. 

We focused on a one-dimensional coronagraph with an un-apodized entrance aperture, 
though other designs may prove useful. For example, the full two-dimensional problem may 
yield surprises. Some coronagraph designs use apodized entrance pupils — though classical 
coronagraphs with apodized entrance pupils generally require M{x) to vary with wavelength 
in a complicated way. It may also be possible to gain benefit by cascading multiple image 
plane masks. But neglecting these avenues, our analysis shows that our existing library of 
one-dimensional image mask designs is complete. 

Thanks to David Spergel, Wes Traub, and Ann Bragg for comments on this manuscript, 
and the folks at the Leiden Coronagraphy Workshop for inspiration. M.J.K. acknowledges 
the support of the Hubble Fellowship Program of the Space Telescope Science Institute. 
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